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We study the effect of the intrinsic (Rashba or Dresselhaus) spin-orbit interaction in 
superconductor-nanowire-superconductor (SNS) weak links in the presence of a spin-splitting field 
that can result either from an intrinsic exchange field or the Zeeman effect of an applied field. 

We solve the full non-linear Usadel equations numericall}f^ and analyze the resulting supercurrent 
through the weak link and the behavior of the density of states in the center of the wire. We point 
out how the presence of the spin-orbit interaction gives rise to a long-range spin triplet supercur¬ 
rent, which remains finite even in the limit of very large exchange fields. In particular, we show how 
rotating the field leads to a sequence of transitions between the 0 and tt states as a function of the 
angle between the exchange field and the spin-orbit field. Simultaneously, the triplet pairing leads 
to a zero-energy peak in the density of states. We proceed by solving the linearized Usadel equa¬ 
tions, showing the correspondence to the solutions of the full equations and detail the emergence 
of the long-range supercurrent components. Our studies are relevant for on-going investigations 
of supercurrent in semiconductor nanowires in the limit of several channels and in the presence of 
disorder. 


I. INTRODUCTION 

The antagonist nature of conventional singlet super¬ 
conductivity and magnetism has been clearly illustrated 
in experiments studying supercnrrents flowing through 
ferromagnetic weak linkP^^. There, the spin-splitting 
(exchange) field h suppresses the supercurrent within 
a typically short magnetic length scale = ^JliD/h, 
where D is the diffusion constant of the wire. As sug¬ 
gested in Ref.m this suppression can be lifted by convert¬ 
ing part of the singlet supercurrent into a triplet with a 
finite projection of the magnetic moment of Cooper pairs, 
by utilizing an inhomogeneous magnetization at the in¬ 
terface between the ferromagnet and the superconduc¬ 
tor. This component couples electrons with spins from 
the same band, and therefore it is not sensitive to the 
spin-splitting field. This suggestion was experimentally 
demonstrated in a number of workJ^ utilizing a series of 
different types of magnetic layers that are non-collinear 
with respect to each other. 

Besides using magnetic materials, the spin-splitting 
field can be realized via the Zeeman effect of an applied 
magnetic field.Similar physics as in the SFS case can 
be envisaged as long as the orbital effect of the mag¬ 
netic field is weak enougtP^ and does not limit the su¬ 
percurrent. Such a situation takes place especially in 
narrow nanowires, where the spin-splitting field in com¬ 
bination with the Rashba-type spin-orbit (SO) interac¬ 
tion has been used in an effort to take these wires to 
the limit of topological superconductivit^fflK^ for the 
detection of Majorana-type excitations at the edges of 
the wires. Most of such experiments are nevertheless 
in the topologically trivial limit. It is hence of interest 
to study the physics of such nanowires in the presence 
of the combination of the spin-orbit and spin-splitting 
fields. This is the aim of the present work. In par¬ 


ticular, we study the supercurrent behavior in systems 
schematically presented in Fig. Contrary to many re¬ 
cent theory works on the effects of spin-orbit coupling on 
proximity superconductivity discussing the fully ballistic 
regimewe assume the wires to be diffusive. Strictly 
speaking this limit requires that all wire dimensions are 
smaller than the elastic mean free path. Typical epitax¬ 
ial nanowires have mean free paths comparable to the 
wire thickness, and much less than the wire length.^^^^^ 
Even in this limit the diffusive-limit theory is likely to 
capture the essential physics much better than the fully 
ballistic limit. In the diffusive limit it is generally possi¬ 
ble to obtain a f ully q uantitative fit with between theory 
and experiments,^^^^ which is why also the quantitative 
details of the theory are relevant. On the other hand, 
the quasiclassical theory we employ corresponds to set¬ 
ting the Fermi wavelength ^ 0. Therefore, it can¬ 
not capture effects related to for exampl e weak a ntilocal- 
ization, possibly relevant in these wiresAlterna¬ 
tive derivations of the quasiclassical theory on the fixed 
few-channel limifP^ cannot be directly connected on the 
many-channel limit considered here. 

This paper extends on the work of Tokatly and 
Bergeretj ^^ l ^^ l who introduced the mechanism of includ¬ 
ing the intrinsic (Rashba or Dresselhaus) spin-orbit in¬ 
teraction as a spin-dependent vector potential into the 
Usadel equation describing inhomogeneous superconduc¬ 
tivity in the diffusive limit. They also pointed out how 
for certain relative orientations of the wire, spin-orbit 
fields and the exchange field, the combination of the lat¬ 
ter two may produce triplet supercurrent that survives 
even at large exchange fields. In particular, they showed 
that the wire with a homogeneous exchange field and in¬ 
trinsic spin-orbit interaction is gauge equivalent to a fer¬ 
romagnet with inhomogeneous magnetization (see also 
Ref. |32j). Here we study this mechanism quantitatively 
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FIG. 1. SNS junction studied in this work: a diffusive 
nanowire of length L connects two bulky superconductors. 
The nanowire is characterized by the spin-orbit field A and 
the system exhibits an exchange field h either due to an ap¬ 
plied magnetic field B or by proximity to magnetic material 
(not in the picture). We assume that both A and h are in 
the xz-plane, and that h is at an angle 0 compared to the (z) 
direction of the wire. 


(Sec. IV). In particular, we show the dependence of the 
supercurrent vs. exchange field for varying magnitudes 
of the Rashba field. We also demonstrate in detail how 
the triplet supercurrent depends on the direction of the 
magnetic field applied in the plane of the wire. For strong 
spin-orbit coupling, we predict that changing the direc¬ 
tion of the field drives the junction through a sequence 
of O-TT transitions. 

Besides supercurrent, we also study the local density 
of states in the junction in Sec.|V| This is also the typical 
observable in the studies of Majorana physics. We find 
that spin-orbit interaction induces a zero-energy peak for 
a range of exchange fields. This peak originates from 
the induced long-range triplet amplitude of the super¬ 
conducting pairing, and it is quite sensitive to the pre¬ 
cise direction of the field and the amplitude of the SO 
coupling. 


II. USADEL EQUATION WITH SPIN-ORBIT 
COUPLING 


to be dominated by the intrinsic SO coupling and neglect 
these other terms in the following. Thus inside the nor¬ 
mal metal we have 

DV • VG^) = [-ie+h -ih-a,G^], (2) 

where = e + ir], and 77 —>■ O'*' is a small term sped- 
fying the location of the poles of the Retarded Green’s 
function. 

Introducing a dimensionless position coordinate z' = 
z/L and defining a Thouless energy Et = where 

L is the length of the normal metal wire, we can work in 
energy units of e = EjE^ and h = \i/Et and use the 
scaled vector potential A^ = AL. Moreover, below we 
assume the Zeeman field to lie in the substrate plane, and 
therefore hy = 0 . Equations 0 and 0 have to be sup¬ 
plemented with the normalization condition (G^)^ = 1. 
In the numerical solutions, we implement this by using 
the Riccati parameterizatiorP^ (see Appendix A) for the 
retarded Green’s function. In that parametrization, the 
Nambu-space Green’s function is specified in terms of 
two parameters, 7 and 7. In general, both of these pa¬ 
rameters are 2x2 matrices in spin space. The Usadel 
equation written for 7 reads 

= 

- 2ie'y - ih ■ [ 7 , a] + [A^ 7 ] + 2{As,j}N{As - ^Ag-y) 
+ 27 (j^dz>j)N{Ag - jAgj) + {Ag - 7 As 7 )V(a^/ 7 )) . 

( 3 ) 

For equilibrium observables, we replace —ie by the Mat- 
subara frequencies ujn = 7rT(2n + l).I^The equation for 
7 is obtained by substituting 7 7, V A and by 

taking a complex conjugate of the scalars. 


We implement the spin-orbit interaction into the Us¬ 
adel equation describing the quasiclassical Nam bn-spin 
retarded Green’s function in the diffusive limifP®SIMI 
(here and below, e = ft = ks = f except when we discuss 
particular values of the observables) 


dv-(g^vg^) 


—isTs — ih ■ c7 -j- A “I" Sgf -1- Sin, G^ 


( 1 ) 

Here D is the diffusion constant inside the nanowire, 5 is 
the energy, a is a vector of Pauli matrices in spin space. 


h = hyTs, h,) and VG^ = VrG^ - i 


Afs, G^ 


IS 


the gauge invariant gradient with A describing the SO 
coupling. The latter is specified below in more detail.l^ 
Inside the normal metal, the superconducting pair po¬ 
tential A is zero and we assume that the term Ein de¬ 
scribing the inelastic scattering is negligible. In practice, 
the junction contains also regular spin-flip scattering and 
scattering due to isotropic spin-orbit coupling,!^ charac¬ 
terized by the self-energy Egf • We assume spin relaxation 


A. Spin-orbit field 

A generic spin-orbit coupling emerging in systems with 
broken inversion symmetry is of the form Ai = aj aj . For 
a one-dimensional wire the symmetry is broken by the ge¬ 
ometry in both directions perpendicular to the wire, and 
therefore the components aj can be considered indepen¬ 
dent of each other. To describe a one-dimensional wire 
in the 2: direction placed on a substrate (see Fig.[^ span¬ 
ning the xz plane, we describe the generic spin-orbit field 
as 


A = {a\ai -h afas, 0 , a^cri -h alaa) . ( 4 ) 

For a thin wire there are gradients only in the z direc¬ 
tion, and therefore the only terms coupling linearly to 
the gradient are related to 0^3. The other terms can 
be non-vanishing and they co ntribu te in general to the 
Dyakonov-Perel spin relaxatiorP^^^, but as they do not 
lead to other interesting physics, we neglect them in our 
numerical results. After neglecting the orbital effect of 
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the field, the properties of the junction depend only on 
the relative direction between h and in • a. 

As we vary the direction of h, we fix Az oc ai, i.e., 
0^3 = 0 . We vary the remaining term a = to ob¬ 
serve the physics due to the intrinsic spin-orbit coupling. 
Note that the term aaiUz gives rise to the linear energy 
term of the form eaaipzlm. Our term a is related to the 
(Rashba) terms of size a considered in the literature on 
Majorana fermions (e.g., 1 ^ by 


where m is the effective mass of the electrons inside 
the nanowire. Our a hence has a dimension of Ti/{eL\ 
where L is some length scale. In the numerics, the im¬ 
portant dimensionless parameter is eaLjh^ where L is 
chosen to be the distance between the two nanowire- 
superconductor contacts. For example with the param¬ 
eters of InSb discussed in Ref. [ 39 ] (m = O.Olbrrie and 
a = 0.2 eV A, where rUe is the electron mass), we get 
ea/Ti ^ A' 10 ^ 1 /m. These estimates are consistent with 
the recent experiments^, which obtain even a somewhat 
larger value of d. For typical wires of length L ^ 1 /im, 
the value of the dimensionless parameter can hence be 
of the order of or larger than unity. On the other hand, 
in order to be able to neglect gradients in transverse di¬ 
rections in Eq. Q, we assume that the wires are narrow 
compared to the spin-orbit length, i.e., wire thickness d 
satisfies d ^ hjea. 


B. Boundary conditions 


C. Supercurrent 

The supercurrent through the junction is characterized 
by the spectral current densit}^ 



( 7 ) 


In the Riccati parameterization, the spectral current den¬ 
sity can be written as 



|Tr 


y'y) N — N (77' — 7^7) N 


{A, 7} jN + A7 {A, 7} N 
+ 7 V {A, 7} jN -h 7V7 {A, 7} TV ] 


( 8 ) 


where j' = d^jdz. 

The supercurrent can be calculated as a weighted av¬ 
erage of the spectral current density 


E / £ \ 

where Rn = L/{Ae^DuF) is the Drude resistance of the 
nanowire in the normal state. Here L is the length, A 
is the area, D the diffusion constant and up the density 
of states at the Fermi level of the nanowire in its nor¬ 
mal state. In the Matsubara technique, the integral can 
be calculated as a sum of the spectral current densities 
evaluated at the Matsubara frequencies ujn as 


For the boundary conditions to the Green’s functions, 
we assume clean NS interfaces. This assumption means 
that the parameters 7 and 7 are continuous across the 
interf ace a nd coincide with bulk BCS superconductor 
value^^^Iini 


i\A\E^ 

iBCS = - / ^2 

5 -h Ty |A|^ - (e + 

IBCS = - , CF2 

£ -h Ty |A|^ - (5 -h T0+)^ 


(6a) 

(6b) 


in the real-time description (for Matsubara frequencies, 
replace e by iuJn)- The exact form of the boundary con¬ 
ditions depends on the chosen form of the Nambu vector. 

In the numerics, we express all the lengths in terms of 
the length L of the wire. In this case, the natural energy 
scale is given by the Thouless energy Et = TiD jI?. We 
mostly concentrate on the limit of long wires, where L ^ 
^0 = \JliD!A. This also means that A/E^ = ^ 

1 . 


Is 


Ef 

cRjsi 


00 

2ttT y] Imjs (uin) ■ 

n=0 


( 10 ) 


In the numerics, we cut the Matsubara sum to the 
index after which the obtained supercurrent changes less 
than 1 %. 


D. Density of states 


The density of states (DOS) is given bjEl 


A^(e,R) = Tr 


G^(e,R)T 3 


( 11 ) 


where Np is the DOS in the absence of superconductivity. 
In the Riccati parameterization, 

N {e, R) = ^ReTr { (1 - 77) + iV (1 - 77)} . (12) 

Besides supercurrent, the density of states is another way 
to characterize the excitation spectrum induced by the 
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proximity effect from the supercurrent. It can be ac¬ 
cessed via a standard tunneling measurement. In Sec. [V| 
we show how a combination of the finite exchange field 
and the Rashba spin-orbit coupling gives rise to a zero- 
energy peak in the density of states. This is qualitatively 
similar to what one expects from the measurements in the 
Majorana wiresalthough the physics of this effect 
is quite different. Initial results for the density of states 
in Rashba wires were presented by us in [42l Recently, 
similar type of results were discussed also in Refs. [43] 
and im but that work concentrated on superconduc¬ 
tor/ferromagnet multilayers in the short junction limit, 
where the junction length is of the order o f the supercon¬ 
ducting coherence length In such mul¬ 

tilayers, the emergence of the long-range triplet super¬ 
conductivity requires either the presence of both Rashba 
and Dresselhaus type spin-orbit coupling, or out-of-plane 
magnetic fields. Moreover, for short junctions the energy 
scales are primarily set by the superconducting gap A in¬ 
stead of the Thouless energy as here. Nevertheless, also 
there the triplet proximity effect leads to the presence of 
zero-energy density of states peaks, and long-range super¬ 
current. However, the quantitative details of the results 
are quite different, and therefore as such not applicable 
for the nanowire setups. 


are given in Appendix B. As a result, we get the spec¬ 
tral supercurrent describing the long-range component. 
In the case of a perpendicular field and (^ = 0), the 
result reduces to 


Im(js) = - 


S'/2a^Dy/u}^ 


(w„ + + |A|2)2 /isinh(\/2VS^) 




(15) 

where 0 is the phase difference between the two super¬ 
conductors, and Cbn = {uJn^‘^Da‘^)lET is the (dimension¬ 
less) Matsubara frequency modified by the pair-breaking 
effect from the spin-orbit coupling i-^ —ie in the real¬ 
time formulation with energy e). Note that the spin-orbit 
coupling a plays here a dual role: first, it induces the 
long-range triplet component of the (spectral) supercur¬ 
rent, and second, it induces pair breaking effects. If we 
would include the ordinary spin-flip or spin-orbit scatter¬ 
ing effects, they would (in the lowest order) induce the 
corresponding terms into Co. 


The observable supercurrent at a given temperature 
T is obtained by summing over the Matsubara fr equen - 
cies or integrated over the real energies as in Eqs. dM - 
At low temperatures ksT % Et we may transform the 
Matsubara sum into an integral and obtain in the long- 
junction limit Et^ e^Tia^D C A (restoring e and ?i) 


III. GENERATION OF THE LONG-RANGE 
TRIPLET COMPONENT 

Below, we present numerical solutions of the supercur¬ 
rent and density of states in a nanowire Josephson weak 
link exhibiting both intrinsic spin-orbit coupling and ex¬ 
change field. To understand these results, let us first 
study what we ex pect to find in the limit of a weak prox¬ 
imity effect which allows us to linearize the Usadel 
equation. 

Linearizing Eq. (§ yields 


Is{T « 0) = 

32ETea‘^D ^_2aLe/h 

hR^Ti 


^ 2aLe 
1 + + 


2a^L^e^ 


sin(0). 

(16) 


Strictly speaking this is valid only for 2eaLlh ^ 1, but 
this approximation fits numerical results reasonably well 
also for low a. On the other hand, at high temperatures 
ksT ^ Et^ it is enough to include only the lowest Mat¬ 
subara frequency and the result is 


din = 2w„7-zh-[7,(7]+[A^,7]+2{As,7}As+2z{a2/7, AJ. 

(13) 

Here we have discarded all the terms 0(7^) and noticed 

that N = 1 due to the normalization condition = 

1. We separate the different spin components as 

3 

7 = (14) 

i=0 

and choose a specific SOC field A^ = aaiUz and direc¬ 
tion of the magnetic field, h = h(0, sin(^), cos(^)). The 
resulting equations describe the interplay between the 
short-range and long-range pairing components. 

We can solve the resulting second-order boundary 
value problem by separating the short-range and long- 
range components and using the previous to find ef¬ 
fective boundary conditions for the latter. The details 


Is{T > Et) = 

__ SErea^D (17) 

{-KkeT + + |A|2)2 hRNh ^ 

where x = 2'iTkBT/E t + 2aLelh. 

We compare the analytical result to the full numerics 
in Eig. H Eor large a ^ 6h/{eL)^ the full theory shows 
a second maximum in the supercurrent, absent in this 
analytical approximation (see Eig. |^. 

Eor a non-zero 0 G]0,7r/2[, both long-range compo¬ 
nents /o and /s become non-zero. Due to their cou¬ 
pling, the supercurrent obtains terms oscillating with 
aLe sm{0)/h. The resulting oscillations vs. 0 for aLe ^ 1 
signal transitions between 0 and tt states, which would be 
observed as cusps in the dependence of the critical cur¬ 
rent on the field direction. This is qualitatively described 
in Appendix B and further explored in Eig. 
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IV. SUPERCURRENT 

The effects of the spin-orbit coupling show up only in 
the presence of a non-zero exchange field in the wire. This 
exchange field can be established either in an intrinsically 
ferromagnetic wire, via the magnetic proximity effect in 
a wire in contact with a ferromagnetic insulator, or via 
the Zeeman field from an applied external magnetic field. 
In all these cases the (direct or stray) magnetic field gen¬ 
erates also the orbital effect for the charge carriers. The 
relevance of t his o rbital effect depends on the aspect ra¬ 
tio of the wir d^ l ^^ l Here we assume that the wire is thin 
enough so that we can disregard the orbital effect and 
concentrate only on the Zeeman field. Unless stated oth¬ 
erwise, the numerical results are obtained for a field in 
the direction of the wire (i.e., 6> = 0). 

The most straightforward experiment is to vary the 
magnitude of the exchange field. In Fig. we plot the 
supercurrent Is at phase difference 0 = 7r/2 (close to 
the critical current) at the temperature T = O.IF^t for 
a long junction (length A/Et = = 1000) as a 

function of the exchange field for varying magnitudes of 
the spin-orbit coupling a. For low a < l/L, the super¬ 
current exhibits damped oscillations as a function of the 
exchange field as shown before for example in Refs.|2]and 
[45l For Is < 0, the junction enters the tt state. With an 
increasing spin-orbit field, the oscillation of the supercur¬ 
rent is supressed and the junction no longer can be found 
in the 7r-state (note that non-quasiclassical cor rectio ns to 
our theory can lead to the presence of a 0-stat^^^Ell with 
0 7 ^ 0,7r due to the spin-orbit coupling). In the absence 
of spin-orbit coupling, a = 0, the supercurrent dies out 
for large exchange fields. However, a finite a yields also 
a finite super current even at rather large exchange fields. 


This minimum persists to rather high temperatures as 
shown in Fig. 



FIG. 3. Supercurrent through the SNS junction as a function 
of the exchange field with different temperatures. 


We can compare our results quantitatively to those ob¬ 
tained from the linearized Usadel equation at large values 
of the exchange field, where the major contribution to the 
supercurrent comes from the long-range components. In 
Fig-i we plot the supercurrent for a few values of the 
exchange field as a function of a. The solid lines show 
the results from the exact numerics, whereas the dashed 
lines come from the Matsubara sum of Eq. (15) (quite 
close to Eq. ®). We find that both approaches yield 
a supercurrent that is non-monotonous wit h respe ct to 
the value of the spin-orbit field. Equations ( 15|16 ) cap¬ 
ture the first oscillation rather well, but the amplitude of 
the exact supercurrent is somewhat larger that obtained 
from the analytics. 



FIG. 2. Exchange field dependence of the supercurrent 
in a SNS junction for different intrinsic spin-orbit coupling 
strengths a. A finite SO coupling strength yields a finite 
supercurrent through the junction even with large exchange 
fields. 



FIG. 4. With large exchange fields, the supercurrent oscillates 
as a function of the spin-orbit field strength reaching a maxi¬ 
mum around a = 2/L. Solid lines show the num erical r esul ts 
and dashed lines are obtained by using Eq. (15) in Eq. (10). 


Eor an intermediate value oi a = 2/L, the supercurrent 
vs. exchange field exhibits a minimum at h ~ 10 ... 20 Lt. 


Tokatly and BergerefP^^^ showed that in order to get a 
long-range triplet supercurrent, the term describing the 
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exchange field should not commute with the vector po¬ 
tential. We check this by investigating the dependence 
of the supercurrent on the angle 0 of an in-plane ex¬ 
change field with respect to the wire direction (see in¬ 
set of Fig. [^. For 6> = 0 the exchange field is in the 
2 ;-direction and therefore produces the term proportional 
to ds, whereas we describe the (Rashba-type) spin-orbit 
term proportional to di. For 0 = 7r/2, both are described 
by terms proportional to di. The supercurrent ls{0) is 
plotted in Fig.[^ Indeed, for 6> = 7r/2, ls{0) only contains 
singlet components of the pairing amplitude and there¬ 
fore almost vanishes because of the large value of the 
chosen exchange field. In addition, as described in Ap¬ 
pendix B, we find that for a ^ l/L, the junction shows a 
sequence of 0 — tt transitions, signalled by the sign change 
of supercurrent at (j) = 7r/2. The magnitude of a dictates 
the position of these transitions. For large exchange field, 
these positions do not depend at all on the field, and they 
only weakly depend on temperature. In experiments, a 
is not usually a variable quantity. However, the angle of 
the applied field can be varied straightforwardly. There¬ 
fore, studying the detailed angular dependence of the su¬ 
percurrent, one is able to determine the magnitude and 
direction of the Rashba vector potential A^. 



FIG. 5. Supercurrent as a function of the angle 0 of the ap¬ 
plied field, h = h(cos 0 a 3 -fsin^ai), h = 170Ft and otherwise 
the same parameters as in Fig. Rotating the exchange 
field makes the long-range component disappear when the 
spin-orbit field is parallel to the exchange field {0 = 7r/2). 
In addition, applying a field at an intermediate angle results 
into an alternating sequence of 0 — tt transitions, their number 
depending on the precise value of a. 


V. DENSITY OF STATES 

Besides supercurrent, the proximity effect from the su¬ 
perconductor on the normal wire can be characterized via 
a tunneling measurement of the density of states. In the 
absence of the exchange field or the spin-orbit field, the 
density of states in the proximity wire exhibits a phase- 


dependent minigap ,1^ whose size can be approximatively 
described via Eg{(j)) ^ 3.1F^t cos^((/)/2). In the absence 
of spin mixing (either via the intrinsic spin-orbit cou¬ 
pling considered here, spin-flip scattering or the isotropic 
spin-orbit coupling), the exchange field simply shifts this 
minigap by ah for spin a = ±. For h > Eg{(j))^ we hence 
expect to see two regions with N{E) = 1/2 in the spin- 
averaged density of states N{E). The intrinsic spin-orbit 
coupling mixes the spins and leads to a closing of the 
minigaps for the spin-resolved density of states. As we 
show below, an intermediate magnitude of the intrinsic 
spin-orbit coupling also leads to a zero-energy peak in 
the density of states. The shape and height of this peak 
is very sensitive to the exact parameters of the system. 
We illustrate these results in the following via a few ex¬ 
amples, a more complete description can be found from 
Ref. |49l 

The density of states is in general position dependent. 
To illustrate the proximity induced effects in a symmetry 
point, we present the results calculated in the middle of 
the normal metal. 



E/Et 


FIG. 6. Local density of states (DOS) in the center of the wire 
with h = S.OOFt and 0 = O.SOtt given with different spin-orbit 
fields. The DOS peaks at zero-energy for a = 2.00/L while 
the peak converts to a zero-energy dip for larger a. 


In Fig. we show the DOS peaking at zero-energy for 
Rashba field oi a = 2lL while for larger a the peak con¬ 
verts to a zero-energy dip and again to an energy gap 
centered around zero energy. As discussed in Refs. l5Ql 
and ESI the particular value for the zero-energy density 
of states results from the competition of the singlet prox¬ 
imity effect aiming to lower ^"(0) and the long-range 
triplet proximity effect to increase it. This results from 
the different symmetry of the singlet vs. triplet compo¬ 
nents. Namely, it is straightforward to show from the 
linearized equations, Eq. ( |B5| ), that for the singlet com¬ 
ponent / 2 (e = 0) = /I (e = 0)5 whereas the triplet compo¬ 
nents satisfy fj{e = 0) = —fj{^ = 0)? / 7 ^ 2. Expanding 
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Eq. (12) to the lowest order in 7 then yields 


N{0) ^ Nf 


1-21/2(0)1' 


2El/i(0)l' 


(18) 


For large /i, the singlet component vanishes on a short 
distance ^ £rn at the interface, whereas the long-range 
triplet components /o ,3 are much larger in the center of 
the wire. As a result, the latter yield A^( 0 ) > 1 . As 
shown in the figures, this may signal the presence of a 
zero-energy peak, but not necessarily. The values of a 
where the transition from a peak to a dip takes place are 
similar to those yielding a maximum in the long-range 
supercurrent. Fig. 




FIG. 7. Density of states with a constant o = 2/L and with 
different values of the exchange field h. The zero-energy peak 
forms only for h ^ 2Eg{(j)). 


Although the origin is different from that expected for 
Majorana junction^^dJti^^ for the diffusive nanowire 
Josephson junctions the zero-energy peak appears only at 
a large enough exchange field, in practice for h ^ 2Eg{(j)) 
and for a restricted range of the values of a. We plot 
the density of states at a constant a for a few different 
exchange fields in Fig.[^ The exact form of the density of 
states depends a lot on the value of the exchange fields, 
and the zero-bias peak decays at large values of h. 

The density of states can also be controlled by applying 
a supercurrent through the junction, so that the phase 
difference (j) across the junction changes. Besides chang¬ 
ing (and closing) the minigap for a vanishing exchange 
field, this changes the form of the density of states. The 
phase dependence is plotted in Fig.|^ Note that typically 
in long junctions (A ^ Et or L ^ = \/^T)/A) the 

density of states for 0 = tt is almost featureless due to the 
destructive interference of the pair amplitudes emanat¬ 
ing from the two superconductors. This is also modified 
by the spin-orbit coupling as shown in Fig.|^ also in this 
case a zero-energy peak forms. 

Let us furthermore demonstrate the connection be¬ 
tween the triplet proximity effect and the zero-energy 


FIG. 8. Phase dependence of the DOS in the middle of the 
nanowire. 



E/Et 


FIG. 9. Even with phase difference of 0 = tt, the DOS peaks 
at zero energy for a finite SO coupling strength a. 


peak.LMlOIIl Therefore, we study the density of states as 
a function of the angle 0 of the exchange field between the 
^ direction of the wire and the x direction, as in Fig. for 
the supercurrent. This is shown in Fig. 10 For 0 = 7 r /2 
we expect to get only the short-range proximity effect. 
In this case we indeed only find the exchange-field split 
minigaps, and the spin-orbit coupling makes almost no 
contribution to the form of the density of states. 


VI. CONCLUSIONS 

In this work we have discussed the detailed effects of 
the Rashba/Dresselhaus -type intrinsic spin-orbit inter¬ 
action on the supercurrent carried through a diffusive 
nanowire. We solve the full Usadel equation to obtain 
the supercurrent of the junction and the local density of 
states. We reproduce the long-range triplet proximity ef¬ 
fect predicted before in the linearized limit and show how 
the resulting supercurrent depends on the direction of the 






































1.4 



FIG. 10. Density of states for different directions of the ex¬ 
change field, specified in terms of the angle 0 (see inset of 

Fig.|5|. 


P+ 


N N-f 




(A2) 

where N = {1 77 ) and TV = (1 + 77 ) . Thus from 


Eq. (Al) the Green’s function is 


= 


N 0 
0 N 


1 — 77 27 

27 -(1 - 77) 


(A3) 


The projectors have the property 


VP± = ±P+ 
where 


VP 


P_ ±P_ 


VP 


+ 5 


(A4) 


applied field. Besides the complicated direction depen¬ 
dence and the predicted zero-energy peaks and dips in the 
local density of states, our results pave the way of quan¬ 
titatively analyzing the experiments carried out in the 
nanowire Josephson junctions. To reach for example the 
topological regime and confirm the Majorana character 
of the excitations in the nanowires, it is important that 
the experimentalists are able to characterize their junc¬ 
tions in detail. Our work gives a fixed point for such a 
characterization, in the (non-topological) many-channel 
diffusive limit. We hence expect this work to be relevant 
as an intermediate step for establishing the experimental 
constraints for using such junctions in topological quan¬ 
tum computing. 
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Timo Hyart for very useful discussions, and Francesco Gi- 
azotto, Attila Geresdi and Charles Marcus for explaining 
the characteristics of InAs/InSb nanowires. This work 
was supported by the Academy of Finland through its 
Center of Excellence program, and by the European Re¬ 
search Council (Grant No. 240362-Heattronics). 


C/=(;;j) and Cf=("o). (A5) 

In the spin-dependent case, 7 and 7 are 2 x 2-spin 

matrices, that is 

3 3 

1 = ^1 3^3 and 7 = ^ 7 ^a^-. (A 6 ) 

i=o i=o 

Using the above relations, we have derived the Usadel 
equation (§ and the spectral supercurrent, Eq. (§. 

Appendix B: Linearized equations 

Linearizing Eq. ^ yields 

7'' = 2cjn7 + '^[7? 7] + 2{As, 7 }As + 2T{7', Ag}, 

(Bl) 

where we use a shorthand notation dz'j = 7 ^ Writing 


Appendix A: Riccati parameterization 


The normalization condition 



1 implies that 


the possible eigenvalues of are ±1. Therefore in spec¬ 
tral representation, G^ can be written in terms of so 
called Shelankov projectors a^^ 


3 

7 = 

i=0 

assuming a general form for the SO coupling 


A = ((aicri + <a2C^35 0? + o^4.o’3) • 


(B2) 


(B3) 


G^ = P+- P- with P± = I (1 ± . (Al) 

The Shelankov projectors and the Green’s func¬ 
tion are convenient to parametrize in the Riccati 
parameterizatiorP^ 


and choosing the exchange field in plane 


h = h{cos0a3 + sinOai) (B4) 

we can separate the equations for the different spin com¬ 
ponents. They read 
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fo = ( 2 w„ + 4 ^ a?)/o + + 4^04/3 (B 5 a) 

/" = (2w„ + 4:{al + al))fi - 2/icos0/2 

+ 4 ( 0 ^ 10^2 + <^3<^4)/3 + 4 ^ 3/0 (B5b) 

f 2 = 2hcos0fi + 2ujnf2 — 2hsm0fs (B5c) 

/g = 4:{aia2 + Q^3<^4)/i + 2hsm0f2 

+ (2cj^ + 4 ( 0^2 + <^4))/3- (B5d) 

and the equations for fi parametrizing 7 are obtained by 
substituting fi ^ fi^ h ^ —h, and by taking a complex 
conjugate of the scalars. 

The boundary conditions read 


/2=ce'^/2j2 = ce-'^/2^c = 


v1^ 


2 ’ 

, 


(B 6 ) 


with (f having the opposite signs at x = 0 and x = L. For 
i = 0,1, 3 the functions vanish at the boundaries, that is 
fi{0) = ML) = Uo) = ML) = 0 . 

Linearizing the spectral current density (see Eq. (§) 
yields 


Im(js) = hm {Tr [(7^ - j'j) - (77' - ^7) (B 7 ) 

+4 {A,7}7 + 7{A,7} +{A,7}7 + 7{A,7} ) 1. 


Using the Pauli matrix expansion from Eq. (B2) and the 
general form for the SO coupling, Eq. (B3) we can sim¬ 
plify Eq. ([B7|) 


lm{M = Im - fjfj) 

-j=0 

+ 4 i((a 3 (/o/i + /i/o) + <^ 4 (/o /3 + /3/0)) 


. (B 8 ) 


In the numerics we have chosen A = (0,0,—adi). On 
the other hand, the analytics becomes more straightfor¬ 
ward by writing the /-function components in the basis 
dictated by the exchange field. Therefore, applying the 
rotation exp(i 6 >cr^/ 2 ) to h, 7 and A yields 


/o =(2w„ + 4 q;^)/o - 4iacos{9)f[ - 4iasm{9)fl^ 

/" =[ 2 w„ + 4 q!^ cos^( 0 )]/i - 2/1/2 - 4 iQ!Cos( 6 >)/o 
+ 2a^ sin( 26»)/3 
M =2hfi + 2uJnf2 

/g =[ 2 w„ + 4a^ sin^( 6>)]/3 + 2a^ sin( 20 )/i - 4 / 0 ; sin( 6 »)/o. 

The corresponding equations for fj are obtained after 
replacing h ^ —h and a ^ —a. 


Eor h ^ a, it is now straightforward to identify fi 
and /2 as the short-range components, decaying within 
the length from the interfaces. These two 

components are the m = 0 triplet and the singlet compo¬ 
nent of the pairing amplitude, respectively. These com¬ 
ponents can thus be solved separately in a straightfor¬ 
ward manner, but the resulting analytic expressions are 
too lengthy to be printed here. 

The short-range components generate boundary con¬ 
ditions for the long-range component s /p and f s that 
decay within the length Iq = 1/^/2ujn + 4(a^), ^3 = 
1 / 1 / 2ujn + 4^2 sin( 6 >) ^ Eirst disregarding the cou¬ 
pling term between /o and / 3 , we get an analytic solution. 


/o = Aq sinh 



fs = As sinh 



+ 4:ia cos{0) / £0 sinh 
Jo 



f[{x)dx 

(B9) 


2a^ sm 


in(2(9) [ 
Jo 


is sinh 


fi{x)dx. 

(BIO) 


This solution takes into account the boundary condi¬ 
tion /o(0) = fs{^) = 0. The prefactors Ao ,3 would 
be obtained from the other boundary condition /o(l) = 
/ 3 ( 0 ) = 1. However, we concentrate only on the vicinity 
of the interface at ^ = 0, and disregard these terms. Note 
that both components vanish for 0 = 7r/2, when the ex¬ 
change field is collinear with the spin-orbit field, whereas 
/3 = 0 for (9 = 0 . 

Eor ^0 ^ ^ ^ Cm? we may expand and get 


h(.z) ? 

^ Aia cos(^) j 

[ fi{x)dx 

0 

(Bll) 

m ^ 

^ 2a^ sin( 2 ^) 

[ fi{x)dx. 

Jo 

(B12) 


The previous expression requires also one partial integra¬ 
tion and using the fact that /i(0) = 0. The equation is 
written for fs{z) because only the derivative tends to a 
constant in this interval. Performing the integral yields 
Jq fi{x) ~ cexp(i 0 / 2 )/( 2 \/h), saturating for ^ ^ Cm- 
Repeating a similar procedure on the other end, z = 1 
(in reduced units) yields finally the full boundary condi¬ 
tions for the long-range components 

/o(0) = 2ia cos{0)ce^^^‘^I\Jh (B13) 

fs{0) = sm{20)ce^^^^/Vh (B14) 

/o(l) = 2iacos{0)ce~^'^^‘^ fVh (B15) 

/^(l) = -0.2 sin(2<9)ce-^^/V\/^, (B16) 


where now z = 0,1 mean the position ^ Cm away from 
the contacts. The boundary conditions for /o ,3 are ob¬ 
tained by changing the sign of the derivatives, a —a, 
(f -H- —(j). The remaining equations for the long-range 
components can be written for the two-component vec- 
tor /= (/o h)^ as 

/" = [ujn + 4q;^]/- 4q!^ cos^(0)(T4,/- 4iasui{9)axf'. 

(B17) 
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We get rid of the second off-diagonal term by defining 
h = U{z)f = ex.p{—2iasm{0)axz)f. This satisfies 

h" = cos^(6>))h — 4(a^ cos‘^{0)U{z)a_iU\z)h 

with h^'\0) = /(O), exp(2i(asin(6>)cra;)h^'^(l) = /(I). Be¬ 


sides Eq. (B18) we can hence transform the coupling 
terms to th e bou ndary conditions. As the general solu¬ 
tion to Eq. (|B18|) is lengthy, we first disregard this term. 


In this case, solving (B18) is straightforward. The full so¬ 
lution for h{z) is lengthy, but the spectral supercurrent 
is given by js = cos^(^)c^ sin(0)/((a) with 


I{a) = 


cos(2(asin(6>)) {2ujn — (a^sin^(6>)) — 4(asin(6>) sin(2(a sm 


^_ 2 e 2 v^\/^ cos(4(a sin(6>)) -h 


(B19) 


where + 2a^cos^{0). In particular, for ^ = 0, 

where the above approximation of neglecting the term 
is not relevant, we get Eq. 

Eor 0 between 0 and 7r/2, the two long-range compo¬ 
nents mix and produce a supercurrent that can change 
sign for a fixed phase as the direction of the magnetic field 
is tuned. Such a sign change is described by the terms of 
the form cos(2(a sin(6>)) that result from the above trans¬ 


formation between / and h. They show that the sign 
change takes place whenever a ^ 1. However, due to the 
approximation of neglecting the term, this result does 
not match very well the exact solution at arbitrary range 
of parameters. Therefore, it is better to use the rather 
straightforward solution of the full linearized equations 
for finding the supercurrent in this case. 


^ The code used for calculating the results in this paper will 
be made available upon the publication of this manuscript 
and before that can be requested from the authors. 
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